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' Abstract 

in ! 

j^. , We find the approximate solutions for the Higgs masses and cou- 

£Nj ' plings in the NMSSM with exact and softly broken PQ-symmetry. 

The obtained solutions indicate that there exists a mass hierarchy in 
the Higgs spectrum which is caused by the stability of the physical 

o 

C^i 1. Introduction 

The minimal SUSY version of the Standard Model (SM) stabilizing the 
mass hierarchy does not provide any explanation for its origin. Indeed the 
Minimal Supersymmetric Standard Model (MSSM) being incorporated in 
the supergravity theories leads to the /U-problem. Within supergravity mod- 
els the full superpotential is usually represented as an expansion in powers 
of observable superfields C a 

W = Wo(h m )+ f i{h m ){H 1 eH 2 ) + h a0J C a C p C J + ... , (1) 

where h m and Wo(h m ) are the "hidden" sector fields and its superpotential 
respectively. The "hidden" sector fields acquire vacuum expectation values 
of the order of Planck scale (Mpi) breaking local supersymmetry and gen- 
erating a set of soft masses and couplings in the observable sector. From 
dimensional considerations one would naturally expect the parameter u to 
be either zero or the Planck scale. If u = then the minimum of the Higgs 
boson potential occurs for (Hi) = and down quarks and charged leptons 
remain massless. In the opposite case, when the values of u ~ Mpi, there 
is no spontaneous breakdown of SU(2) x U(l) symmetry at all since the 
Higgs scalars get a huge positive contribution u 2 to their squared masses. 
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In order to provide the correct pattern of electroweak symmetry breaking, 
\i is required to be of the order of the electroweak scale. 

In the simplest extension of the MSSM, the Next-to— Minimal Super- 
symmetric Standard Model (NMSSM) ^H], the superpotential is invariant 
with respect to the discrete transformations C' a = e 2m / 3 C a of the Z 3 group. 
The term fj,(HiH 2 ) does not meet this requirement. Therefore it is replaced 
in the superpotential by 

W H = XSiHteHz) + , (2) 

where S is an additional superfield which is a singlet with respect to SU(2) 
and U(l) gauge transformations. A spontaneous breakdown of the elec- 
troweak symmetry leads to the emergence of the vacuum expectation value 
of the extra singlet field (S) = sj\pl and an effective ^-term is generated 
(// = Xs/y2). The Z% symmetry of the superpotential naturally arises in 
string inspired models, where all observable fields are massless in the limit 
of unbroken supersymmetry. 

In this paper we investigate the Higgs masses and their couplings to 
the Z-boson in the NMSSM using approximate solutions. In Section 2 we 
specify the Higgs sector of the model. In section 3 the exact Peccei-Quinn 
(PQ) symmetry limit in the NMSSM is studied and approximate solutions 
for the Higgs masses and couplings are obtained. The scenario of soft PQ- 
symmetry breaking is discussed in section 4. In section 5 we summarize 
our results . 



2. NMSSM Higgs sector 



The NMSSM Higgs sector involves two Higgs doublets H\£ and one singlet 
field S. The interactions of the extra complex scalar S with other particles 
is defined by the superpotential (J2j) that leads to a Higgs boson potential 
of the following form: 



V 



,/2 



l^ll 2 



|2\2 



+X 2 \S\ 2 {\H 1 \ 2 + \H 2 \ 2 ) + Xk 



S* 2 {H l eH 2 ) + h.c. 



+ k 2 |5| 4 + 



(3) 



XA x S{H x eH 2 ) + -A K S 3 + h.c. 



+X 2 \(H 1 eH 2 )\ 2 

+m 2 l \H 1 \ 2 + m 2 \H 2 \ 2 + m 2 \S\ 2 + AV 



+ 



where g and g' are SU(2) and U(l) gauge couplings respectively, while AV 
corresponds to the contribution of loop corrections. The couplings g, g' , A 



2 



and k do not violate supersymmetry. The set of soft SUSY breaking pa- 
rameters includes soft masses m 2 , m|, nig and trilinear couplings A K , A\. 

At the physical minimum of the potential © the neutral components of 
the Higgs doublets H\ and H2 develop vacuum expectation values V\ and 
V2 breaking the electroweak symmetry down to U(l). Upon the breakdown 
of SU(2) x U(l) symmetry three goldstone modes (G^ and G°) emerge, 
and are absorbed by the and Z bosons. In the field space basis rotated 
by an angle (3 (tan/3 = V2/V1) with respect to the initial direction 

Im H$ = (P sin + G° cos 0)/y/2, H{ = G~ cos (3 + H~ sin /3 , 
ImH$ = (P cos (3-G° sin ff)/y/2, H+ = H+ cos (5 - G+ sin (5 , (4) 
ImS = P s /V2 

these unphysical degrees of freedom are removed by a gauge transformation 
and the mass terms in the Higgs boson potential can be written as follows 



V m ass = M 2 H± H + H- + 1(P P S )M 2 ( * j+±(Hh N)M 2 \ h 



H 

h 

N 



(5) 

where we replace the real parts of the neutral components of the Higgs 
doublets by their superpositions H , h so that 

ReH^ = (hcos (3 - H sin(3 + V!)/V2 , ReS = (s + N)/V2, 

(6) 

ReH$ = (hsin (3 + H cos (3 + v 2 )/V2 . 

^ , P , fdV dV dV \ „ , 

From the conditions for the extrema — — = — — = — — = of the 

\OVl OV2 OS ) 

Higgs effective potential © one can express m/a, mf, m 2 . via other funda- 
mental parameters, tg/3 and s. Substituting the obtained relations for the 
soft masses in the 2x2 CP-odd mass matrix M 2 - we get: 

T2 __2 _ V (_ « 



M 22 = -^-3 + Sm 2/9 ~ 3 T^ + A 22 , ( 7 ) 

2 2 A 

M, 2 , = Ml = V2Xvii ( — X — — 2— ) + A 



12 = ^21 = V^av/z -— - - + z\ M 
smzp A/ 

As 1 / K 



where u = V^T+^f = 246 GeV, /i = —=, x = — [A\ + 2-//) sin 2/3 and 

V2 2/i V A / 

Ajj are contributions of the loop corrections to the mass matrix elements. 

The mass matrix (J7J) can be easily diagonalized via a rotation of the fields 

P and P s by an angle 6 A (tan26» j4 = 2M 1 2 2 /(M 1 2 1 - M| 2 )) . 
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The charged Higgs fields H ^ are already physical mass eigenstates with 

\2 V 2 

M 2 H± =m 2 A - — + M 2 v + A ± . (8) 

Here My/ = is the charged W-boson mass and A± includes loop cor- 
rections to the charged Higgs masses. 

In the rotated basis H , h , N the matrix elements of the 3x3 mass 
matrix of the CP-even Higgs sector can be written as P]— @]: 



I 

A 2 



Ml x = m\+ (j-y J « 2 sin 3 20 + An 



M| 2 = Mf cos 2 2(3 + — v 2 sin 2 2(3 + A 22 

r2 _ X„2 , K A „ , A ^ 2 kA „ ,2 



M 33 = 4^ Z + -A K /i+— X-— ^Sill2 i a + A33, (9) 

M 2 2 = M| 1= ^-^^ 2 sin4/? + A 12 , 

Mf 3 = M| 1 = - V^Av^x cot 2/3 + A i3 , 

M| 3 = M| 2 = y / 2Au/i(l - x) + A 23 , 



where Mz = —v is the Z-boson mass, g = g 2 + g' 2 , and Ay can be 

calculated by differentiating AV 3 . Since the minimal eigenvalue of a 
matrix does not exceed its smallest diagonal element, at least one Higgs 
scalar in the CP-even sector has to be comparatively light: < y M 22 . 
At the tree level the upper bound on the lightest Higgs mass in the NMSSM 
was found in 0. It differs from the corresponding theoretical limit in the 
minimal SUSY model only for moderate values of tan (3. As in the MSSM 
the loop corrections from the t-quark and its superpartners raise the value 
of the upper bound on the lightest Higgs mass resulting in a rather strict 
restriction on m/^ < 135 GeV j^j. The Higgs sector of the NMSSM and 
loop corrections to it were studied in jSJ. 

In the field basis P, Ps, H, h, N the trilinear part of the Lagrangian, 
which is responsible for the interaction of the neutral Higgs states with the 
Z-boson, is simplified: 



Lazh = ^MzZ^Z^h + -Z/j, 



H(d„P) - (d^H)P 



(10) 



Only one CP-even component, h, couples to a pair of Z-bosons while an- 
other, H, interacts with pseudoscalar P and Z. The coupling of h to the 
Z pair is exactly the same as in the SM. In the Yukawa interactions with 
fermions h also manifests itself as the SM like Higgs boson. 
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The couplings of the physical Higgs scalars to the Z pair (gzzi, i = 
1, 2, 3) and to the Higgs pseudoscalars and Z boson {gzA\i and 9zA 2 i) appear 
due to the mixing of h, H and P with other components of the CP-odd 
and CP-even Higgs sectors. Following the traditional notation we define 

the normalized .R-couplings as: gzzi = Rzzi x |Mz and g Z Ap = TjRzAji- 
All relative -R-couplings vary from zero to unity and are given by 



Rzzi = > RzAn = -Ujji sin 9 A , RzA 2 i = cos 6 a , (H) 



where is unitary matrix relating components of the field basis H, h, and 
N to the physical CP-even Higgs eigenstates. 

3. Exact Peccei— Quinn symmetry limit 

First of all let us discuss the NMSSM with k = 0. At the tree level 
the Higgs masses and couplings in this model depend on four parameters: 
A, fj,, tan (3, m^or x). When A is small enough (say A < 0.1) the experi- 
mental constraints on the SUSY parameters obtained in the minimal SUSY 
model remain valid in the the NMSSM. If tan/3 < 2.5 the predominant part 
of the NMSSM parameter space is excluded by unsuccessful Higgs searches. 
Non-observation of charginos at LEPII restricts the effective //-term from 
below: \pb\ > 90 — 100 GeV. Combining these limits one gets a useful lower 
bound on wia at the tree level: 



Requirement of the validity of perturbation theory up to the high energy 
scales constrains the parameter space further. In order to prevent the 
appearance of Landau pole during the evolution of the Yukawa couplings 
from the electroweak scale to Grand Unification scale (Mx) the value of A 
has to be always smaller than 0.7 . 

In the NMSSM with k = the mass of the lightest pseudoscalar van- 
ishes. This is a manifestation of the enlarged SU(2) x [C/(l)] 2 global sym- 
metry of the Lagrangian. The extra U(l) (Peccei-Quinn) symmetry is 
spontaneously broken giving rise to a massless Goldstone boson (axion) [Jj. 
The Peccei-Quinn symmetry and the axion allow one to avoid the strong 
CP problem, eliminating the #-term in QCD 8 . At low energies the axion 
gains a small mass due to mixing with the pion. The mass of orthogonal 
superposition of P and Ps is 



m\ > 9M z x . 



(12) 



2 2 i 

m A2 =m A + 



2 



X + A; 



22 • 



(13) 
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The lower bound on tua (|12|) leads to the hierarchical structure of the 
CP-even Higgs mass matrix. It can be written as 

eCt \ 

2 ' (14) 
e 2 B J 



M 




where e << 1. Indeed the top-left entry (M^ 1 = A) of the corresponding 
3x3 mass matrix © is the largest one in the dominant part of parameter 
space. It is proportional to m 2 A while M 2 
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M| 2 



Ml 



M§ and 



M 2 3 ~ ra^Mz- Therefore the ratio Mz/rriA plays the role of a small 
parameter e. 

The CP-even Higgs mass matrix can be reduced to block diagonal form: 



VM 2 V ] 



M ii + 



Mf 3 



M 2 X 



\ 



0(e 3 ) 
0(s 3 ) 



0(e 3 



M 2 2 



0(e 3 



M 



23 



M 



23 



M? 3 M 2 2 
M 2 , 



M 



33 



by virtue of unitary transformation 



V 



( 



\ 



£ 



2 

-er 



r t r 



2 



-rrt 



M 2 U 
M 2 U 



CA~ 



(15) 



(16) 



The obtained matrix (|15|) is now easily diagonalized via a rotation by an 
angle 9 of the two lowest states 



R 




tg29 



M| 3 



M? 3 M 2 2 
M 2 X 



M| 2 



Mj 3 + 



M 2 n 



(17) 



As a result we find the approximate formulae for the masses of the CP-even 
Higgs bosons 



m 



h 3 



h 2 ,hi 



,,2 M 13 



M| 2 + M| 3 




Mfi 



M 13 

33 M x 2 i 



+ 4 



*4 



M? 3 M 2 2 
M 2 U 



(18) 
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Also using the explicit form of the unitary matrix IP ~ R', that links 
H, h and N to the mass eigenstates, the approximate expressions for the 
couplings of the lightest Higgs particles to the Z-boson can be established: 



Rzz2 ~ cos 8 



Rzzi ~ - sin 8 



Rza x 2 



RzA Y \ 




sm 



sm 



8 A 



8a 



(19) 



The approximate solutions for the CP-even Higgs boson masses and 
couplings shed light on their behaviour as the NMSSM parameters are var- 
ied. As evident from Eq, (fl8|) . at large values of tan (3 or fi the mass-squared 
of the lightest Higgs scalar tends to be negative because M| 3 becomes large 
while bottom-right entry of the matrix (|15|) goes to zero. Due to the vac- 
uum stability requirement, which implies the positivity of the mass-squared 
of all Higgs particles, the auxiliary variable x is localized near unity. At 
the tree level we get 



where m A = 2/u/sin2/3. The allowed range of the auxiliary variable x is 
quite narrow (see [3]). According to the definition of tua the tight 
bounds on x enforce tua to be confined in the vicinity of ^tg(3 which is 
considerably larger than the Z-boson mass. As a result the masses of the 
charged Higgs boson, heaviest CP-odd and CP-even Higgs states are rather 
close to mA- At the tree level the theoretical bounds on the masses of the 
lightest Higgs scalars are 



The masses of hi and h\ are set by the Z-boson mass and Xv respectively 
so that m/j 1 , mh 2 <C tua in the allowed range of the parameter space. 





rn 



(20) 



'A 




(21) 
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4. Soft breaking of the PQ-symmetry 



Unfortunately searches for massless pseudoscalar and light scalar particles 
exclude any choice of the parameters in the NMSSM with k = 0, unless 
one allows A to become very small (U]. In order to get a reliable pattern 
for the Higgs masses and couplings the Peccei-Quinn symmetry must be 
broken. Recently different origins of extra U(l) symmetry breaking were 
discussed |1U|-|11|. Here we assume that the violation of the Peccei-Quinn 
symmetry is caused by non-zero value of k. As follows from the explicit 
form of the mass matrices J7J), (JSJ) and Q in this case the Higgs spectrum 
depends on six parameters at the tree level: A, k, //, tgf3, A K and wia (or x). 
We restrict our consideration by small values of n when the PQ-symmetry 
is only slightly broken. To be precise we consider such values of n that do 
not change much the vacuum energy density, the last requirement places a 
strong bound on k when A goes to zero: 

k < A 2 . (22) 

If k 3> A 2 then the terms 1 5*| 4 and — A K S 3 in the Higgs effective potential 

3 

(j3J) becomes much larger \fj,\ 4 ~ M| increasing the absolute value of the 
vacuum energy density significantly. A small ratio n/X may naturally arise 
from the renormalization group flow of A and k from Mx to Mz @J ■ 

The soft breaking of the PQ-symmetry does not lead to the realignment 
of the Higgs spectrum preserving its mass hierarchy. Still M 2 X is the largest 
matrix element of the CP-even Higgs mass matrix in the admissible part of 
the NMSSM parameter space. Therefore the approximate formulae (|18|) - 
(|19j) obtained in the previous section remain valid in the considered limit. 
It is easy to see that the lightest CP-even Higgs states respect a sum rule 

m 2 hl + m 2 h2 = Ml + Ml - ^| . (23) 



The right-hand side of Eq. (|2.3|) is almost insensitive to the choice of mA and 
rather weakly varies with changing tg(3. As a result of the sum rule (|23|) 
the second lightest Higgs scalar mass is maximized as m/ ll goes to zero, and 
vice versa the lightest Higgs scalar gets maximal mass when m/^ attains 
minimum (see also Fig.l). According to Eq. (j!8|) the mass of the lightest 
CP-even Higgs boson vary within the limits: 

< ml < mm|M 2 2 ,M 2 3 -^| (24) 

The mass matrix of the CP-odd Higgs sector also exhibits the hierar- 
chical structure. Indeed the entry M\ x is determined by m\ whereas the 
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off-diagonal element of the matrix (J7J) is of the order of Xv ■ tua- Since 
the ratio k/X is small, the other diagonal entry M22 m \- this again 
permits one to seek the eigenvalues of the matrix Q as an expansion in 
powers of Xv/tua- The perturbation theory being applied for its diagonal- 
ization results in concise expressions for the squared masses of the Higgs 
pseudoscalars 



M 2 +^ 
M22 + Ml, 



Ml 



Ml 



(25) 



Because the PQ-symmetry is now broken the lightest CP-odd Higgs boson 
also gains non-zero mass. In compliance with Eq. ()19|) the couplings of h\ 
and /12 to a Z-pair obey the following sum rule: 



R-ZZl + R-ZZ2 — 1 ' 



(26) 



while Rza x \ and Rza-ci are suppressed by a factor (X 2 v 2 /m A ). 

At the tree level and large values of tg/3 (tgfl > 10) the approximate 
expressions (|18j) and (|25jl describing the masses of the Higgs scalar and 
pseudoscalar particles are simplified 



m 



h 3 



2 2 A 2 U 2 

m A 2 = m A + —2~ x ■ 



-3-A k /j,, 



in 



h-2, hi 



M\ + 4^/j 2 + jA Kf i± 



(27) 



±\ M 2 -4^ 2 -^) +8A2t,V(l-x) 2 



A 2 



A 



making their analysis more transparent. Again the positivity of the mass- 
squared of the lightest Higgs scalar restricts the allowed range of x 



V2kM z 



X 2 i 



< x < 1 + 



s/2kM z 



X 2 v 



(28) 



if the PQ-symmetry is only slightly broken. When k<A 2 the admissible 
interval of the auxiliary variable x shrinks drastically establishing a very 
stringent bound on the value of tjia and strong correlation between mA, H 
and tg (3 (see ^^) that constrains the masses of the heavy Higgs bosons 
m/^ m H ± rriA 2 in the vicinity of \i tgf3. 

The results of the numerical studies of the Higgs boson masses and their 
couplings including leading one-loop corrections from the top and stop 
loops are given in Fig. 1-3. As a representative example we fix the Yukawa 
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couplings at the Grand Unification scale so that X(Mx) = k(Mx) = 
2h t (M x ) = 1.6, that corresponds to tg(3 > 3, \(M t ) ~ 0.6 and n{M t ) ~ 
0.36 at the electroweak scale. We set fx = 150 GeV which is quite close to 
the current limit on n in the MSSM. The parameter A K occurs in the right- 
hand side of the sum rule (|23|) and in the mass of the lightest pseudoscalar 
m Ai with opposite sign. As a consequence whereas m 2 ^ rises over changing 
A K , m 2 Al diminishes and vice versa. Too large positive and negative values 
of A K pull the mass-squared of either lightest scalar or pseudoscalar below 
zero destabilizing the vacuum that restricts A K from below and above. To 
represent the results of our numerical analysis the parameter A K is taken 
to be near the center of the admissible interval A R ~ 135 GeV. 




200 300 400 500 600 700 



m A (GeV) 

Figure l.The dependence of the neutral Higgs boson masses on uia for 
\ = 0.6, k = 0.36, n = 150 Ge V, tgf3 = 3 and A K = 135 GeV. The 
one-loop Higgs masses of scalars (solid curve) and pseudoscalars (dashed- 
dotted curve) are confronted with the approximate solutions (dashed and 
dotted curves). Masses are in GeV. 

In Fig. 1-3 the masses of the neutral Higgs particles and their couplings 
to Z are examined as a function of mA- From the restrictions 1)28(1 on the 
parameter x and numerical results presented in Fig.l it is evident that the 
requirement of the stability of the physical vacuum and the experimental 
constraints on ji and tg(5 rules out low values of m^, maintaining mass 
hierarchy whilst k < A 2 . The lightest Higgs scalar and pseudoscalar can be 
heavy enough to escape their production at LEP. Moreover as one can see 
from Fig. 2 the lightest Higgs scalar can be predominantly a singlet field, 
making its detection more difficult than in the SM or MSSM. The lightest 
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Higgs pseudoscalar is also singlet dominated, making its observation at fu- 
ture colliders quite problematic; the coupling of the lightest CP-even Higgs 
boson to a CP-odd Higgs bosons and a Z is always strongly suppressed (see 
Fig. 3) according to Q19|) . The hierarchical structure of the mass matrices 
ensures that the heaviest CP-even and CP-odd Higgs bosons are predom- 
inantly composed of H and P. As a result the coupling Rza 2 3 1S rather 
close to unity while Rzz3 is almost negligible. In Fig. 1-3 the approximate 
solutions (|18 [) - ()19j) are also given. They work remarkably well. 




200 300 400 500 600 700 



m A (GeV) 

Figure 2. The absolute values of Rzzi and Rzz2, plotted as a function of 
niA for the same values of X, k, fi, tan/3 and A K as in Fig.l. Solid and 
dashed-dotted curves reproduce the dependence of Rzzi and Rzz2 on rriA 
while dashed and dotted curves represent their approximate solutions. 



5. Conclusions 

In the present article we have obtained the approximate solutions for the 
Higgs masses and couplings in the NMSSM with exact and softly broken 
PQ-symmetry which describe the numerical solutions with high accuracy. 
The approximate formulae (|18|) - (|19|) provide nice insight into mass hier- 
archies in the considered model. The vacuum stability requirements and 
LEP restrictions on the NMSSM parameters leads to the splitting in the 
spectrum of the Higgs bosons. When k = or k < A 2 the charged Higgs 
states, the heaviest scalar and pseudoscalar are nearly degenerate around 
m,A ~ fJ-tdft- The masses of new scalar and pseudoscalar states, which 
are predominantly singlet fields, are governed by the combination of pa- 
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rameters — fi. In the NMSSM with exact and softly broken PQ-symmetry 
A 

they are considerably lighter than the heaviest Higgs states. Decreasing 
k pushes their masses down so that they can be even the lightest parti- 
cles in the Higgs boson spectrum. The SM like Higgs boson has a mass 
around 130 GeV. We have established useful sum rules for the masses of 
the lightest Higgs scalars and their couplings to a Z pair. Also we found 
that the couplings of the lightest CP-even Higgs states to the lightest pseu- 
doscalar and Z-boson are suppressed. Observing two light scalar and one 
pseudoscalar Higgs particles but no charged Higgs boson, at future colliders 
would present an opportunity to distinguish the NMSSM with softly broken 
PQ-symmetry from the MSSM even if the heavy states are inaccessible. 




200 300 400 500 600 700 

m A (GeV) 

Figure 3. The absolute values ofRzA x \ and Rza x 2 as a function ofrriA- The 
parameters X, n, fx, tan (3 and A K are taken to be the same as in Fig.l. Solid 
and dashed-dotted curves correspond to Rza x \ and Rza x 2 while dashed and 
dotted curves represent their approximate solutions. 
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